clear
use Figures02_03.dta, replace

twoway ///
	(scatter RESP_L DATE, mcolor(black) msize(vsmall) msymbol(o)) ///
	(scatter RESP_C DATE, mcolor(black) msize(vsmall) msymbol(t)) ///
	(scatter CONT_L DATE, mcolor(gs5) msize(vsmall) msymbol(o)) ///
	(scatter CONT_C DATE, mcolor(gs5) msize(vsmall) msymbol(t)) ///
	(scatter COOP_L DATE, mcolor(white) msize(vsmall) msymbol(o)) ///
	(scatter COOP_C DATE, mcolor(white) msize(vsmall) msymbol(t)) ///
	(lowess response DATE, lcolor(black) lwidth(medthick) lpattern(solid)) ///
	(lowess contact DATE, lcolor(gs5) lwidth(medthick) lpattern(solid)) ///
	(lowess cooperation DATE, lcolor(white) lwidth(medthick) lpattern(solid)) ///
	(scatter RESP_L DATE if DATE==0, mcolor(black) msize(medium) msymbol(O) connect(l) lcolor(black) lwidth(medthick) lpattern(solid)) ///
	(scatter RESP_C DATE if DATE==0, mcolor(black) msize(medium) msymbol(T) connect(l) lcolor(black) lwidth(medthick) lpattern(solid)) ///
	(scatter CONT_L DATE if DATE==0, mcolor(gs5) msize(medium) msymbol(O) connect(l) lcolor(gs5) lwidth(medthick) lpattern(solid)) ///
	(scatter CONT_C DATE if DATE==0, mcolor(gs5) msize(medium) msymbol(T) connect(l) lcolor(gs5) lwidth(medthick) lpattern(solid)) ///
	(scatter COOP_L DATE if DATE==0, mcolor(white) msize(medium) msymbol(O) connect(l) lcolor(white) lwidth(medthick) lpattern(solid)) ///
	(scatter COOP_C DATE if DATE==0, mcolor(white) msize(medium) msymbol(T) connect(l) lcolor(white) lwidth(medthick) lpattern(solid)) ///
	, ///
	legend( order(50 "Response:" 10 "Landline" 11 "Cell" 51 "Contact:" 12 "Landline" 13 "Cell" 52 "Cooperation:" 14 "Landline" 15 "Cell") ///
		rows(3) region(lcolor(none) fcolor(gs10))) ///
	plotregion(fcolor(gs10)) ///
	ytitle(Rate) ylabel(0(20)100, angle(0) nogrid) ///
	xtitle("") xlabel(16071 16802 17532 18263 18993 19724 20454 21185,angle(0) labsize(vsmall)) ///
	scheme(s1mono) 
graph export Figure03.pdf, replace
